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Abstract 

We construct biorthogonal polynomials for a measure over the complex 
plane which consists in the exponential of a potential —V {z, z") and in a 
set of external sources at the numerator and at the denominator. We use 
the pseudonorm of these polynomials to calculate the resolvent integral 
for correlation functions of traces of powers of complex matrices (under 
certain conditions). 



1 Introduction 

In several domains of physics like, for instance, in string theory [1] and 
in quantum chromodynamics [2] we are interested in the spectral properties of 
random matrices; most results in this domain are included in the correlation 
functions for the invariants of the matrices hke determinants, traces... Although 
most interesting results are expected from infinitely large matrices, exact results 
can be obtained in appropriate cases for finite matrices using the technique of 
orthogonal polynomials. The appropriate cases are when the integrals over the 
matrix elements can be transformed into integrals over the eigenvalues of the 
matrices; this is the case for hermitian matrices where a unitary transformation 
reduces the integrals to the real line for each eigenvalue. When the matrices 
are complex, there exists a small class of potentials [3] (which includes the 
Gaussian potential) where the integration can be reduced to the complex plane 



(or a compact domain of the complex plane) for each eigenvalue. Moreover this 
is valid only for a restricted class of correlation functions which is essentially 
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in a potential V (M, M+) , which becomes 
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in a potential V {z, z*) . In this context, given a potential V {z, z*) such that the 
partition function is defined as 
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we wish to calculate the following integrals 
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In the special case where Li = Mi and L2 = M2 the application of 
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n (~af^) n (~'^) ^ defines the resolvent 
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Clearly enough, the large ^j, 77* expansion of is a formal power series the 
coefficients of which are the correlation functions (2) . 

In this publication, we calculate 7jv and Jn; the formal power series of Jn 
will be done somewhere else. We now describe the method which is based on the 
existence of biorthogonal polynomials for the potential V {z, z*) in the presence 
of external sources. 

We consider a real potential V {z, z*) which admits an infinite set {pn {z)} 
of orthogonal monic polynomials 



J (fz p*^ {Z) Pn {Z) e-^(^'^*) = hn Sr, 



(6) 
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where (z) is a short notation for [pm (z)]* . More generaUy, we introduce a 
positive Borel measure /x {z, z*) on the complex plane and write (6) as 

J dll {Z, Z*) {z) Pn {z) = hn 6nm (7) 

We introduce four operations: 

operationl : d^^ (z, z*) = (z — C) {z, z*) 

operation2 : dfi2 {z, z*) = {z* — rf) djj {z, z*) 

operations : d^^ {z,z*) = d^{z,z*) 

operation^ : dji^ {z, z*) = —rz^ dji {z, z*) 

Although these four measures are not real, each of them admits an infinite 

set of monic biorthogonal polynomials. The successive iterations of the four 
operations construct an infinite set of monic biorthogonal polynomials for the 
measure 

dji {z, z*\ £,i, r]*;yi, x*) = g d^ {z, z*) (8) 

f[{z*-x*) f[{z-yi) 

i=l i=l 

These monic biorthogonal polynomials are denoted by qn{z',£,i,'r]*;yi,x*) and 
In (2;; Vi, d; Xi, y*) and satisfy 

/..(^.^•;«.,,.-;».,.:) C(.;,.,S;..,».-) ,„(.;{.,,.•;».,<) Hl,»f 

(9) 

Clearly, g* {z; rj^, Ci'i^i^ V*) is deduced from g„ (z; ^i,r]*;yi, x*) by complex con- 
jugation and by the exchange rj^^ yi <==^ Xi which implies Li <;=^ 
L2, Ml ^=4> M2. The unicity of the biorthogonal polynomials requires that we 
stay away from the surfaces of the (^j, r]*;yi, x*) space where \\qn\\ = 0. 

This construction is a generalization of the well-known result by Christoffel 
[4] which shows that, given a positive Borel measure of one variable dfi (x) on the 
real line and its infinite set of orthogonal polynomials, it is possible to construct 
an infinite set of othogonal polynomials for the measure 

L 

dti{x;Q=l[{x-^,) dfi{x) (10) 

i=l 

The integrals In corresponding to this measure have been calculated by Brezin 
and Hikami [5] . 
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Then, Uvarov [6] in 69 and recently Fyodorov and Strahov [7 — 8] constructed 
an infinite set of orthogonal polynomials for the measure 



i=l 



M 



dfi {x) 



(11) 
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if the variables yi are chosen away from the real axis. In (2003), Akemann 
and Vernizzi [9] (see also ref.[10]) constructed an infinite set of (bi)orthogonal 
polynomials corresponding to the measure 



L2 



Li 



dfx{z,z*;^-,Vi) = ]l(z*-rj*) ]l{z-^i) dii{z,z*) 



(12) 



i=i 



where dji {z, z*) is a positive Borel measure over the complex plane. 

In the following sections we construct the biorthogonal polynomials corre- 
sponding to the measure^ (8). We now explain how to calculate the integrals In 
and Jjv from the biorthogonal polynomials. As usual, we write the expression 



A (2) = n - 



%<3 



TTJV-I (2^1) ••• T^N-l [zn) 
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(13) 



where the polynomials 7r„ (z) are any set of monic polynomials. If moreover, 
the polynomials 7r„ {z) arc (bi)orthogonal, the variables of integration factorize 
in the integrals (3 — 4) and the (bi) orthogonality can be used to obtain 



Af-l 
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(14) 

^! n (15) 

i=0 



i=0 



In the general case, the pseudonorms are given in section 3; in the special 

i > (16) 



case where Li = Mi and L2 = M2, we have 

II Il2 I. 



A-1 

where £), is defined in (19) and D-i is given in (108-109); consequently 

Dn-1 



IN 



(17) 
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or from (108-109) 



M2 
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n i^t-rj*) n im-^j) 

Ml(Mi-l) M2(Af2-l) j_ j-l i j — l 

= ^"^ ' ^"^ ' A (a.*) A(y) A(r) A(0 ^^-^ 



(18) 



where A (?/) is the Vandermonde determinant of the variables yi and £)„ is a 
determinant expressed in terms of four kernels Kn, Nn, An, N* defined in 
section 2 



Dn 



Nn{S,l,yMi) ... N,i{(.MiiyMi) An{xl,yMi) 



■ An{x*M,,yi) 

■ N*{xm2,Vm2) 

(19) 



Although A>, (S,,y) has a single pole at ^ = y, this pole is cancelled by the 
prefactor in Ipf. Let us mention that this type of determinant has also been 
found in Ref. [11] for the correlation functions of several hermitian matrices 
coupled by a chain of potentials. Finally, the resolvant Jjv is calculated in 
section 4 as the determinant 



J. 



where "r>„"is 



"det" 



N 



JV-l 



I^n (Cmi : '?Mi ) An(r]\, ) 
Kn{£.M,,Vl) Hn iVl,Vl) 

Kn (^Mi > »?M2 ) l*?! , Vm^ ) 



(20) 



-^n ('7m2 ' ^Mi ) 
II n (^M2 5 ^Ma) 

(21) 



The kernel Hn ($, y) is defined in (26) and concists in subtracting the pole 
of Nn (C, y) at ^ = y. The determinant "-D„" is obtained from the determinant 
Dn by changing all 7V„ on the diagonal into Hn, changing all y^ into ^j., all Xk 
into Tj^ and finally ignoring all double poles at = and at 77* = 77^ as we 
develop the determinant (which we denote by " "). In fact, there is no single 
poles either at = ^j. or at ?7| = r]f. since the residues are zero. 

We finally mention that the large ^j, r/^ expansion of Jn (which will be 

done somewhere else [12]) is related, in the case of the Gaussian potential, to 
the number of graphs obtained from (1) by Wick's contraction. The large N 
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expansion of this number classifies the graphs according to their genus (t'Hooft 
expansion ); another large N limit which is considered in [12] is the so called 
BMN limit [13 - 14 - 15] where, in (1), -j^ and are kept constant. 



2 The functions and the kernels 



In this section wc describe the biorthogonal polynomials for the four operations 
described in the introduction; then we give several simple examples of measures 
before describing the general case and the proofs in section 3. 



Pn (z) 
P*n{z) 
K 



operationl: 



g„(2:;^, $;$,$) = 



1 



g:(^;$,r;$,^)= K{z;e) 

Pn+l (C) 



Pn+l {z) Pn+l iO 
Pn {z) Pn (C) 

hn 



Pn (0 



PniO 



Win] 



-hn 



PniO 



(22) 



In (22), the symbol $ means an empty set of variables and the kernel {z; ^*) 
is 

Pt (z) Pi (0 



Kn{z,n = J2- 



i=0 



hi 



(23) 



The iteration of operationl to several sources is the generalization to po- 
tentials of two variables of the Christoffel's construction with potential of one 
variable (in the case of one real variable Qn {x; ^) and g'* {x; ^) have similar 
forms as in (22) and are equal by Christoffel-Darboux relation; such a relation 
does not exist with two variables z and z*). 



Pn {z) 

Pn (z) 
hn 



- operation2: 

qn{z;^,r]*;^,^) = Kn{z;r)*) 
1 



ri* 



Pn iv) 

P*n+l (z) Pn+l iv) 
Pn (Z) Pn iv) 



Pniv) 



\Qn\\ 



-hn 



Pn+l M 
Pn iv) 



(24) 
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- operations: 
We introduce the function 



tn (2/) = j 



djj. {z, z*) pI {z) 



z-y 



(25) 



which, in the case of potentials of one variable is called the Cauchy-Hilbert 
transform of the polynomial Pn {x) . We also introduce the kernels 



Hn {z,y) 
Nn{z,y) 



E 

1=0 

1 



y- z 



Pi{z) tjjy) 
hi 

+ Hn {z, y) 



(26) 
(27) 



Then, 
Pn{z) 
Pn (z) 



q„{z;^,^;y,^) = {z-y) Nn-i{z,y) 



tn-i (y) 



n > 



gf;(2;;$,$;$,2/* 



Pn (z) tn {y) 

Pn-l (Z) tn-1 (y) 



I ||2 , tn (y) „ 

UnW =-hn-i- n>0, 

tn-1 [y) 



190 I 



tn-1 {y) 
= to (y) 



n>0 
(28) 



Pn {z) 
Pn (z) 
hn 



operation4: 



Pn (z) t*n (x) 

Pn-1 (z) t*n_^ {x) 



qn{z;^,^;^,x*) = 

q*n{z-^,^-X,^) = {z* -X*) N*n_^{z,x) 



||g„f = -/i„_i n>0, \\qo\ 



t*n-l {X) 



tn-1 i^) 

hn-1 
t*n-l (X) 

= to (x) 



n > 
n > 



(29) 



We see from these four operations that the biorthogonal polynomials take 

different forms according to the measure. Before proceeding to the general case 
and to the general proofs, let us give some simple examples which make the 
general result easier to understand. We introduce the function 



Q{x*,y)= [ dn{z,z*) — 1— 

J [z* — x*) {z — 



y) 



and the kernel 



(30) 



(31) 



i=0 
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We consider the measure 

dii{z,z*\£,,rf\y,x* 
Then, it is easy to find that 



{z* -X*) {z-y) 



d/i (z, z*) 



(32) 



Pn {z) 



Pn (z) 



z-y 

z-i 



D„ = 



Pn {0 tn i^) 

Nn{z,y) Nn{^,y) An{x*,y) 

Kn{z,v*) K„{^,V*) N*ir],x) 
N^{^,y) Anix*,y) 
Kn{^,V*) Kirj,x) 



Dn-l 



n > 



= K 



n > 0, 



|9o| 



ho {x* - V*) {y - Do (33) 



We let the reader calculate g* (z; rj, £,*;x, y*) from g„ (z; r]*;y, x*) by complex 
conjugation and by the exchange ^ <^ r], x y. 

Now, if we let the variable r;* ^ cx) in (32-33), we obtain the biorthogonal 
polynomials for the measure 



dii{z,z*;C,<^;y,x*) 



{z-0 



djj {z, z*) 



[z* — X*) [z — y) 
Performing the corresponding expansion in the kernels, we obtain 



(34) 



Pn {z) 



D„ = 



qniz;^,<^;y,x*) 

Pn (z) Pn (?) t;, (x) 

Pn-1 (z) Pn-1 (?) (x) 

N„iz,y) N^i^y) A„{x*,y) 

Pn (?) C (x) 

Nn{^,y) Ar^{x*,y) 



Z-y 
z-^ 



Dn- 



n > 



hn\ 



-hn-1 



Dn 
D. 



n-1 



n>0, llgof = -(y-0 i^o 



(35) 



If we let the variable ? ^ oo in (32-33), we obtain the biorthogonal polyno- 
mials for the measure 

dii {z, 0*; ?, y, x*) = -[^^^^ d^^ {z, z*) (36) 

(z* - X*) {z - y) 

Performing the corresponding expansion in the kernels, we obtain 



8 



Pn {z) 



= (z-y) 



Dr, = 



\Qn\\ 



Qn {z;^,'^>;y,x*) 

Nn-i{z,y) An-i{x*,y) 

tn{y) An{x*,y) 



hn-1 
Dn-1 



n>0 



-K- 



D. 



n-1 



n>0, \\qo\\' = {x* - T,*) Do (37) 



Also, if we let in (32-33) the variable a;* — > oo and if we use the following 
result proved in Appendix A 



tn (y) 
An {x*,y) 



K 

— as w ^ cx) 

yU+l ^ 

Pn+l (0 

as V ^ 00 

yn+2 » 

c+1 (^) 

— , „ as V — > c)0 

yn+2 » 



we obtain the biorthogonal polynomials for the measure 

dn{z,z*;^,7i*;y,<i>)^iz* -rj*) ^-^—^ d^{z,z*) 



{z-y) 



(38) 
(39) 
(40) 

(41) 



as 



Pn{z) qn{z;^,ri*;y,^) 



Nniz,y) Nni^,y) 

Kn{z,TJ*) Kn{^,V*) 



= -hr, 



Z-y 
z-^ 

tn+l{y) Nn{^,y) 
Pn+liv) Kn{^,V*) 

D 



Dn-l 



n > 



Dr. 



n>0, \\qo\\ =-ho {y - ^ £>o (42) 



Finally, if we let in (32-33) the variable y — > cxd we obtain the biorthogonal 
polynomials for the measure 



dfi{z,z*;^,r]*;y,^) = 



I'l dn{z,z*) 
[z* — x*j 



(43) 



as 



Pn{z) qn{z;^,T]*;^,x* 



z-i 



Pn+l (z) Pn+l {Cl C+1 {x) 
Pn {z) Pn (0 t*n {x) 

Kniz,V*) Kn{^,V*) N*iTJ,x) 



Dn- 



n > 
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Dr,. = 



= h 



1 — n > 0, 



n > 



\\qof = -ho {x*-r,*) Do 



(44) 



The rest of this section is devoted to six other simple examples: 



1°) 



which gives 



(45) 



Pn{z) q„{z;^,^;y,^) 



z-y 



Pn {z) Pn (0 

iV„(z,2/) Nr^itv) 



1 



Nn^iity) 



n > 



Nnii,y) 
Nn-i{^,y) 



n>0, \\qo\\=ho No{^,y) (46) 



2°) 



which gives 



rfM {z, z*; V*; X*) = S^* ''2 df^ 



{z* — X*) 



(47) 



Pn{z) -> g„(0;C,$;y,$) 



z-y 

z-i 



Pn {Z) tl {X) 

K„{z,rj*) N*iv,x) 



1 



N:-iiv,x) 



n > 



h„. n>0, \\qof=ho {x*-r,*) iV*(r?,x) (48) 



3°) 



which gives 



K-i{v,x) 

dn {z, z*; ^, 77*; $) = {z* - v*) {z - dn {z, z*) 



(49) 



Pn{z) g„(0;C,77*;$,$) 



z-i 



Pn+l (z) Pn+l (0 



1 



n > 



4°) 



II l|2 u Kn+1 {^,r}*) 

1 

dfj.{z,z*;^,^;y,x*) 



{z* — X*) {z — y) 



dji {z, z*) 



(50) 



(51) 
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which gives 

Dn 



Pn-l {z) (x) 

Nn-i{z,y) An-i{x*,y) 



(z-y) 



An-i{x*,y) n > 0, Dq^ -Q{x* ,y) 



Dr. 



Dn-1 



n> 0, 



llaoir 



-Dn 



which gives 

Pn (z) 



dii {z, z*; X*) = ^\ dii {z, z*) 



{z* — X*) 



qn{z;^,^;<l>,x*) 



1 



Dr, = 



Pn+1 (Z) Pn+1 iO tn+1 (x) 

Pn (z) Pn iO Ki {x) 

Pn-l {z) Pn-1 iO t*^-l {x) 

Pn+1 iO tn+1 {x) 
Pn (0 tn (X) 



Dn 
Dn-1 



n>0, 



Dr 



-Dn 



6°) 

which gives 

Pn {z) 



dfi {z, z*; i>, T]*;y, $) = ^ ^ dn {z, z*) 

{z-y) 



(z-y) 



Dn = 



qn{z;^,^;^,x*) 

Nn-l{z,y) tniy) hn-1 

Kn-i{z,r]*) pl{ri) Dn-1 

tn {y) tn+1 {y) 

Pniv) Pn+liv) 

-hn-1 n > 0, \\qof = -Dq 

Dn-1 



n>0 
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3 The biorthogonal polynomials 



We now consider the measure (8). We learned from the previous section how to 
construct the biorthogonal polynomials; clearly, it is of the form 



Ml 

Qn [z; £,i,Vi; Vi, Xi ) = — 

ri(^-^i) 



(57) 



where is a factor which makes the polynomial Qn monic. The form of the 
determinant depends of two parameters 

d = n + Li-Mi (58) 
7 = n + L2-M2 (59) 

and its construction is always related to the linear properties of its columns. 

Ml 

Clearly, the expression JJ (-2 — Z/i) * 1. 1^ 7 ™ust be a polynomial in z of degree 

i=l 

(n + Li) which vanishes when z = ^^ for any i. In the determinant, the first 
column is the only z dependant column; then, some of the Li columns at the 
right of the first one are chosen to be identical to the first one with z successively 

Ml 

replaced by ^i, Consequently, the polynomial (2 — Vi)* \-\d,-y vanishes 

i=l 

for z = for any i, and can be divided by JJ (^; — to give a polynomial of 

i=l 

degree n. 

If > we may introduce in the first column the polynomial pd (z) , the 
kernels Kd {z, rj*) , Nd {z, y) ■ li d < it is convenient to introduce 

Nk<o{z,y) = — (60) 

z-y 

Kk<o{z,V*) = (61) 

Ak<oix*,y) = -Q{x*,y) (62) 

Ml 

so that if = — 1 the expression JJ (2^ — -^-1 {z, yk) is a polynomial in z of 

i=l 

degree (n + Li) . Finally, if < — 1 we introduce in the determinant a line 

I N_i{z,y) po{y) ... p_d-2 (y) \ z^ as z ^ 00 (63) 
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The remaining columns are constructed in such a way that 

dij{z,z*;£,i,T]*;yi,x*) pI^{z) qn{z;£,i,r]*;yt,x*) =0 for m < n 
or equivalently 



i=i 

M2 



i=l 



-viz,z') ^0 for m < 



The square pseudonorm is calculated by 
j (fzp*^{z) g 

n(^*-<) 

i=l 



We use the property 



i=l 

M2 

n(^*-<) 



where 



and where the polynomial 



2* —X* 

L2 



i=l 



m+L2-M2 

7rm+i2-M2 (^*) = ^ ai p* {z) 

1=0 



Pm(^M2)yM2(''i) 
^*~^M2 
PM2-2 (2!M2) 

Po (a;M2) 



Olm+L2-M2 = 1 



is monic of degree {m + L2 — M2) in z* and is if (m + L2 — -^2) < 0. 
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If the first column of the determinant contains a term ph (z) , the integrals 
(65-66) generate a term 



Pk {z) 



I J 



+ 



A(a;*) 

P*m iV*i)tl (Xl) 

Po (xi) 



P*miXM2)fM2 iVi)tl {XM2 
P*M2-2 {XM2) 

Po{xm2) 



(70) 

where ak hk requires < fc < (m + £2 — M2) and otherwise. Similarly, the 
integrals (65-66) generate the terms 



m+L2-M2 



Nk{z,y) 



A(a;*) 



i=SMp(0,fc+l) 

P*m{Xl)fl{ri*)Ak{xl,y) ... P*m{XM2)VM2il*i) 

P*M2-2iXl) ... P*M2-2i^M2) 



Kk{z,'ni) 



Po (.Xl) 
m+L2-M2 

i=Sup{a,k+i) 



PoixM2) 



A(x*) 



P*mixi)'Pl{vl)N^{r]i,Xi,) ... P*m{XM2)'PM2ivt)N^iVl,XM2) 
P*M2-2ixi) ... P*M2-2{XM2) 



Po (xi) 



Po{xm2) 



where the sum J2 is absent ii k > {m + L2 — M2) or if (m + L2 — ^2) < 0. The 
equations (71-72) are still valid if < 0. 



3.1 7 = n + L2-M2<0 

Since, in that case, there is no polynomial nm+L2-M2 (z*) in (67) for m <n, the 
calculation of (65-66) as given by (70-71-72), may be written 



Pk (z) 
Nk{z,y) 
Kk{z,Vk) 



t*k (x,) 
Ak {x*,y) 
Nk {Vk^Xi) 



(73) 



where the coefhcients /?j are the same in the three cases. Consequently, a de 
terminant containing 

Pk {z) t*k {Xi) 

Nk{z,y) Ak{x*,y) (74) 
Kk{z,V*k) Nl{r)k,Xi) 
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gives zero in (65) and (66) for m < n. Now, we observe that the coefficients /3j 
are such that 



V,^ ^P*ci^^)\_^ ifc<-(m + L2-M2 + l) , 
^MPc(^0J"11 ifc=-(m + L2-M2 + l) I 

i=l 



a°) if (i > 0, then <;„ {z',^i,r]*;yi,x*) contains a determinant of size 
(ii + M2 + 1) X (is + Ml + - 7 + 1) 



9n {z]^i,Vi;yi,x*) 



Ml 
i=l 



Pd {z) Pd i^j) 



Po{z) Poi^j) t^{xi) 

Op* (xi) 

{xi) 

Nd{z,yi) Nd{^J,y^) Ad{x*,y,) 

Kd{z,vl) Kd{^^,vl) N^iVk,xi) 



Dd,, 



(76) 



where Dd^-y is written in Appendix B ((Bl) for > and (B2) for d = 0). The 
pseudonorm is generated in (76) from the terms (a;;) ; from (75) we find 



Dd, 



(77) 



where -Dd+1,0 is given in (B4). 



b°) if d = —1, then g„ {z; ^j, ry|; yi, x^) contains a determinant of size 
(ii + M2 + 1) X (i2 + Ml - 7) 



(ln{z;£,i,r)*;yux*) 



X{{z-y,) 

i=l 



-1 











p%_^ {xi) 

N_i{z,yi) N_i{^.,yi) A.^ix^yi) 
Nl,{r]k,xi) 



(78) 



where -D-1,7 is given in (B3). The pseudonorm is stU obtained from the tej:mp'^^_^ 
as 

(79) 



1,7 
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where -Do,7+i is given in (B2) if 7 < — 1 and in (B5) if 7 = —1. 



c°) if d < —1, then Qn {z;^f,r}*;yi,x*) contains a determinant of size 
(Li +M2-d)x (L2 + Ml - 7) 

Ml 



i=l 













N-i{z,yi) poiUi) 




pS (xi) 

(a;,) 

p-d-2{yi) N_i{^j,yi) A_i{x*,y,) 

iV*i(r?fe,a:0 



(80) 



where Dd^-y is given in (B3) The pseudo norm is still obtained from the term 
P--y-i (xi) as 



Dd+i, 



7+1 



(81) 



where -Dd+1,0 is given in (B8). 



3.2 7 = n + L2-M2> 

a°) we suppose > 7. 

In this case, for m < n, equations (65-66) give from (70-71-72) 

Pk{z) ] M, ( tl{Xi) \ 

Nd{z,y) }^J2f^^{ ^'ii^hy) } fc>7 + l (82) 

Kd{z,vl) J -1 [ N^{Vk,Xi) J 

while these equations give for = 7 

".<"} (83) 

[ K + Ei=i /3i (a^i) m = n \ 
Consequently, the polynomial is a determinant of size 
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(ii + Ms + 1) X (i2 + Ml + d - 7 + 1) 

Pd (z) 

11 [z - Vi) 



1=1 



Pi {z) Pi (Cj) (a;;) 
Nd {z, Vi) Nd (Cj , Vz) Ad {xj, Vi) 
Kd{z,vl) Kd{^^,vl) N^iVk,^i) 



Ddn 



(84) 



where Dd^-y is given in (B4) if d > 7. and in (B5) if d = 7. The pseudonorm 
comes from the term Pj (z) and is found from (83) to be 



II l|2 / \d-t , -D(i+1,7+1 
hnW = (-) /i7 

where Dd+1,^+1 is also given in (B4) if > 7. and in (B5) if d = 7. 
b°) we suppose = 7 — 1 
In this case, for m < n, we still have in (65) 



Nd{z,y) 
Kd{z,V*k) 



Ah 



i=l 



Ad {x*,y) 
^d iVk^x^) 



however, for m = n 



Nd (z, v) 
Kd{z,vl) 



p;{v) 



Consequently, the polynomials contain a determinant of size 

(Li +M2 + 1) X (L2 + Mi) 

9n {z;£,i,v*i;yi,x*) 

Ml 

^ i=i Nd {z, yi) Nd , Vi) Ad {x^ , ?;,) 

Kd{z,r,l) Kd{^j,vl) N*iVk.xi) 

i=l 

for d > -1 

(ln{z;£,i,ri*;yi,x*) 

7V_i(^,y,) A^-i(^,-,yO A_^{x\,yi) 
A^*i(r?fc,arO 



Dd,i 



(85) 



(86) 



1^' { } 



(88) 



for d 



Y^^z-u) 

i=l 
-1 



(-1) 
£>-l,0 

(89) 
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where is given in (B6) if d > —1 and -D_i_o is given in (B8). The 

pseudonorm is obtained from (87) as 



Ik" 



- hd 
Do,i 



for d 



d> -1 



-1 



(90) 
(91) 



-D-1,0 

where D^+i^^+i is also given in (B6). 

c°) we suppose — l<(i<7 — 1 ^7>0 

We are now in the domain where equations (65-66) give from (70-71-72) 



Nd{z,y) 



m— n+7 

E 

i=d+l 



ti {y) 

P*i iVk) 



Ms 
i=l 



Ad {x*,y) 
NdiVk^Xi) 



(92) 



This relation leads to biorthogonal polynomials which contains a determinant 
of size 

{L-i + M2+j-d)x {L2 + Ml) 



Qn {z;£,i,-q*;yi,x*) 



Ml 

1=1 \.~) 



7-d-l 



hd 



Li 



i=l 



Nd{z,yi) i[7-i,d+i] (yj) Ad{x'^,yt) 
Kd{z,V*k) Pl^-i,d+iiiVk) Kd{ipr]l) N*{r]^,xi) 



for d > 



(93) 



Qn {z;^i,Vi;yi,Xi) = 



Ml 

i=l 

i=l 



(-) 



7-1 



-D-1,7 



for d 
where 



N_i{z,yi) t[7-i,o](yi) N_i{^j,yi) A^i{x*,y^) 
P[7_i,o] ('Jfc) N*_,{rjk,xi) 

-1 



t[7_i,d+i] (y,;) = t^-i (?/,:) t~,-2{yi) ■■■td+i{yi) 

P[-y-l,d+l]{Vk) = P*t-l{Vk) Pj-2{Vk)-Pd+liVk) 



(94) 



(95) 
(96) 
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and Dd^-y is given in (B6) for d > —1 and -D-1,7 is given in (B7). The pseudonorm 
is found from (92), with to = n, as 

Iknf = (-r-^/i.^^ d>-l (97) 

Iknf = ^ d=-l (98) 

where D^+i^^+i is also given in (B6). 



d°) we suppose d < —1, 7 > 



Equations (92) are stiU valid but we have to apply (63). The determinant 
is of size 

(Li + M2 + 7 - d) X {L2 + Ml) 



Ml 



9n {z;^i,r]*;yi,x*) 



N-i{z,yi) t[^-i.o]{yi) P[o~d-2]{yi) N^i{Cpyi) A_l{x*,y^) 
P[7-i,ol(%) iV*i(%,a;() 



for 7 > 



(99) 



Ml 

l[(^-yi) 



Qn iz;^i,r]*;yi,x*) 



(-)" 



7V_i(2;,yi) p[o-d-2](yi) iV-i (^^-,^4) .4_i(x;,?/,) 
iVli(r7;^,a;;) 



for 7 = 



where 



(100) 



P[o -d-2] (yi) = Po ivi) Pi (yi) •••P-d-2 ivi) (101) 

and where D^,-/ is given in (B7) if 7 > and D^.o is given in (B8). Finally, the 
pseudonorm in that case is obtained from (92), with m = n, as 



(102) 
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where Dd+i^^+i is given in (B7). 



We note that the pseudonorm can be written, for all values of d, and 7 as 



qnf = V7 (103) 



with the convention 



hd,i = hinf{dn) d 7 > 

hd^j = — 1 d or 7 < 

hd,-y = 1 d and 7 < (104) 



4 The integrals 1^ and Jn 



As shown in (14-15), the expression for 1^ is related to the product of the 
pseudonorms from to TV — 1. Since the pseudonorms are essentially the ratio 
of two determinants , we obtain from (58-59) 

N-l 

Y[ hi + Li-Mi.i+ L2-M2 
_ ^_^N(Li-Mi-L2+M2) i=0 



D 



N+Li-AIi, A+L2-M2 



N-l 



i=0 



D 



Li-Mi,L2-M2 



(105) 

with the convention (104) and where the determinants £'li-Mi,L2-M2 and 
Z?]v+Li-Mi, N+L2-M2 are given in the previous section depending of the val- 
ues of A^", Li — Ml and L2 — -M2. In the special case where Li = Mi and 
L2 = M2 the integrals In are simply 



D 



N-l 



IN 



(106) 



where 



Dn-1 
D-i 



= Dn,n 



I) 



0,0 



Nn-1 {^j,yi) An-1 {xi,yi) 

KN-l{^j,V*k) N*^-l{Vk^Xl) 

7V_i(e,,2/^) A_i(x*,y,) 
A*i(77fe,Xi) 



iV_i(e,-,yi) iV^i (77^,3;;) 



(107) 



(108) 
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with 



1 

Vl-il 

1 



1 

1 

Vmi — 



A(y) AjQ 



The result (106-107-108-109) is nothing but equation (18). 



In order to obtain the resolvent Jjv from 7jv, we must perform the derivatives 

M2 Ml 



M2 

n 



_d_ 



Ml 

n 



Wi) A{x*) A{y) A(r?*) A(0 



(110) 

This operation can be decomposed into the derivatives in involving the 
Ml first column of Dn while the derivatives in ij* involve the M2 last rows of 
Dn- Both operations are similar and we describe here the derivatives in We 
define 



A(t/) A(0 



(111) 



while the factors {yi — are distributed over the columns of Given I c 
we write 



n(-A)«(,„,)]^^^jn(-|-){n<»-«.)-»); 



yi=ii 
(112) 

The right [.] is easily calculated since for the columns i ^ I we obtain 1 on the 
diagonal and otherwise, and for the columns i € I we obtain Hn (Ci,^,) on 
the diagonal, the remainder being imchangcd except for yk = Cfc- Wc call Z?„ (J) 
the corresponding subdeterminant with the indices i € I. Now, we calculate the 
left [.] ; for any J C (^1, ■■■,^Mi) have 



= if card (J) is odd 



(113) 



and if card(J) is even 



(-)" 



Ml (Ml— 1) 



E n 

all pairings in J {j,k) 



(114) 
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Altogether, we have obtained 



(-)" 



Ml (Ml ~1) 



E 

I even 



E n 

all pairings in J 



-1 



Dn (/) 



Mi(Mi-l) 

(-) ' "A, 



(115) 



where is the determinant obtained from _D„ after the following changes: 

the diagonal is Hn (^j, ^j) , all variables yi arc changed into ^j, finally " " means 
that all double poles at = arc ignored as wc develop the determinant. In 
fact, there is no single poles at ^j. = either, since the residues are zero. 



The same result applies for the derivatives (~sfr) ^^'^ gives 
Jn as in (20). 



5 Appendix A 



We compute the large y behaviour of the function t„ (y) and of the kernels 
Nn (^, y) and An {x*,y) defined in (27) and in (31). The formal power series for 
tn {y) when y ^ oo can be written as 



tn{y) = -E ^ (Ai) 

p=0 ^ 

dfiiz,z*) P*n{z) zP (A2) 

Clearly, Oj, = for p < n. Let us write 

p 

2^ = ^ cCpA; pfe (z) , apk = if fc > p, = 1 (A3) 



Consequently, 



tn {y) = -h„ "Pn (A4) 



tn iv) ~ as 2/ ^ cx) (A5) 
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Now, if we insert the expansion (A4) for f„ (y) and if we use (A3), we obtain 

^ Pn (0 tn jy) _ 1 , . 

71=0 " s y 

If this expression is reported in the definition of Nn (^, y) we obtain 

NAU) = - E (A7) 

p=n+l P 

^r.{^,y) ~ ^^^^ as 2/ ^00 (A8) 
Similarly, we find for Q {x* ,y) defined in (30) 

Q(.%y) = E^^4^ (A9) 

so that 

.^2, t* (x) tr, (y) 
Ar,{x*,y) = - E ' \ (AlO) 

p— n+1 ^ 

~ as 2/ ^00 (All) 



6 Appendix B 

We collect the various forms of the determinant D^^^ according to the values of 
d and 7. 



6.1 7<0 



la°) d > 



Kd- 









to i^l) 

Po i^i) 

P-^-1 (xi) 

Ad-1 {x*i,yi) 

Nd-iiVk,xi) 



(Bl) 
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lb°) d = 



Do, 



lc°) d<0 






Po iVi) 








Po {xi) 



P^-i (xi) 

Nl,{r]k,xi) 











Po (xi) 



(a;;) 

P-d-iiVi) iV-i(^j-,2/i) A_i{x^,yi) 
A^*i(%,arO 



6.2 7 > 

2a°) d>7 



-Dd,7 = 



Pd-i (e,) 



2b°) d = 7 
2c°) < d < 7 



Kd-i{^j,V*k) ]^d-iiVk,Xi) 
Nd-i{ij,yi) Ad_i{x*i,yi) 



t[7-i,d] (yi) Nd (Cj , yt) Ad {x*i , yi) 



where {yi) andpj'^^_^j (jj^) are defined in (95-96). 



2d°) rf< -1, 7 > 1 



^[7-1,0] (y*) P[o-d-i]{yi) N^i{i^,yi) A_i{x*i,yi) 
p\^-^,o]{r}k) iV*i(r/fc,x;) 

where P[_d_i,o] (j/i) is defined in (101). 

2e°) d<-l, 7 = 

„ _ P[o-d-i]{yi) N_i{^j,yi) A^i{x*,y,) 
N*_,{Vk,xi) 
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